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Current OFDMA systems group resource blocks into subband to form the basic feedback unit. 
Homogeneous feedback design with a common subband size is not aware of the heterogeneous channel 
statistics among users. Under a general correlated channel model, we demonstrate the gain of matching the 
subband size to the underlying channel statistics motivating heterogeneous feedback design with different 
subband sizes and feedback resources across clusters of users. Employing the best-M partial feedback 
^ ' strategy, users with smaller subband size would convey more partial feedback to match the frequency 

' selectivity. In order to develop an analytical framework to investigate the impact of partial feedback 

o : 

Xf^ , and potential imperfections, we leverage the multi-cluster subband fading model. The perfect feedback 

scenario is thoroughly analyzed, and the closed form expression for the average sum rate is derived for 
CO ' the heterogeneous partial feedback system. We proceed to examine the effect of imperfections due to 

channel estimation error and feedback delay, which leads to additional consideration of system outage. 
. . , Two transmission strategies: the fix rate and the variable rate, are considered for the outage analysis. 

X ■ 

^ . We also investigate how to adapt to the imperfections in order to maximize the average goodput under 
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heterogeneous partial feedback. 
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I. Introduction 

Leveraging feedback to obtain the channel state information at the transmitter (CSIT) enables a wireless 
system to adapt its transmission strategy to the varying wireless environment. The growing number of 
wireless users, as well as their increasing demands for higher data rate services impose a significant burden 
on the feedback link. In particular for OFDMA systems, which have emerged as the core technology 
in 4G and future wireless systems, full CSIT feedback may become prohibitive because of the large 
number of resource blocks. This motivates more efficient feedback design approaches in order to achieve 
performance comparable to a full CSIT system with reduced feedback. In the recent years, considerable 
work and effort has been focused on limited or partial feedback design, e.g., see HI and the references 
therein. To the best of our knowledge, most of the existing partial feedback designs are homogeneous, 
i.e., users' feedback consumptions do not adapt to the underlying channel statistics. In this paper, we 
propose and analyze a heterogeneous feedback design, which aligns users' feedback needs to the statistical 
properties of their wireless environments. 

Current homogeneous feedback design in OFDMA systems groups the resource blocks into subband 
El which forms the basic scheduling and feedback unit. Since the subband granularity is determined by 
the frequency selectivity, or the coherence bandwidth of the underlying channel, it would be beneficial to 
adjust the subband size of different users according to their channel statistics. Empirical measurements and 
analysis from the channel modeling field have shown that the root mean square (RMS) delay spread which 
is closely related to the coherence bandwidth, is both location and environment dependent Q, 0. The 
typical RMS delay spread for an indoor environment in WLAN does not exceed hundreds of nanoseconds; 
whereas in the outdoor environment of a cellular system, it can be up to several microseconds. Intuitively, 
users with lower RMS delay spread could model their channel with a larger subband size and require less 
feedback resource than the users with higher RMS delay spread. Herein, we investigate this heterogeneous 
feedback design in a multiuser opportunistic scheduling framework where the system favors the user 
with the best channel condition to exploit multiuser diversity Q, lH. There are two major existing 
partial feedback strategies for opportunistic scheduling, one is based on thresholding where each user 
provides one bit of feedback per subband to indicate whether or not the particular channel gain exceeds 
a predetermined or optimized threshold ifTI- llTOII . The other promising strategy currently considered in 
practical systems such as LTE [11] is the best-M strategy, where the receivers order and convey the M 
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best channels llT2l - |[T9l . The best-M partial feedback strategy is embedded in the proposed heterogeneous 
feedback framework. Apart from the requirement of partial feedback to save feedback resource, the study 
of imperfections is also important to understand the effect of channel estimation error and feedback delay 
on the heterogeneous feedback framework. These imperfections are also considered in our work. 

A. Focus and Contributions of the Paper 

An important step towards heterogeneous feedback design is leveraging the "match" among coherence 
bandwidth, subband size and partial feedback. Under a given amount of partial feedback, if the subband 
size is much larger than the coherence bandwidth, then multiple independent channels could exist within 
a subband and the subband-based feedback could only be a coarse representative of the channels. On the 
other hand, if the subband size is much smaller than the coherence bandwidth, then channels in adjacent 
subbands are likely to be highly correlated and requiring feedback on adjacent subbands could be a 
waste of resource; or a small amount of subband-based partial feedback may not be enough to reflect 
the channel quality. In order to support this heterogeneous framework, we first consider the scenario of 
a general correlated channel model with one cluster of users with the same coherence bandwidth. The 
subband size is adjustable and each user employs the best-M partial feedback strategy to convey the M 
best channel quality information (CQI) which is defined to be the subband average rate. The simulation 
result shows that a suitable chosen subband size yields higher average sum rate under partial feedback 
conforming the aforementioned intuition. This motivates the design of heterogeneous feedback to "match" 
the subband size to the coherence bandwidth. The above-mentioned study, though closely reflects the 
relevant mechanism, is not analytically tractable due to two main reasons. Firstly, the general correlated 
channel model complicates the statistical analysis of the CQI. Secondly, the use of subband average rate 
as CQI makes it difficult to analyze the multi-cluster scenario. Therefore, a simplified generic channel 
model is needed that balances the competing needs of analytical tractability and practical relevance. 

In order to facilitate analysis, a subband fading channel model is developed that generahzes the 
widely used frequency domain block fading channel model. The subband fading model is suited for 
the multi-cluster analysis. According to the subband fading model, the channel frequency selectivity is 
flat within each subband, and independent across subbands. Since the subband sizes are different across 
different clusters, the number of independent channels are heterogeneous across clusters and this yields 
heterogeneous partial feedback design. Another benefit of the subband fading model is that the CQI 
becomes the channel gain and thus facilitate further statistical analysis. Under the multi-cluster subband 
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fading mode|3 and the assumption of perfect feedback, we derive a closed form expression for the average 
sum rate. Additionally, we approximate the sum rate ratio for heterogeneous design, i.e., the ratio of the 
average sum rate obtained by a partial feedback scheme to that achieved by a full feedback scheme, in 
order to choose different best-M for users with different coherence bandwidth. We also compare and 
demonstrate the potential of the proposed heterogeneous feedback design against the homogeneous case 
under the same feedback constraint in our simulation study. 

The average sum rate helps in understanding the system performance with perfect feedback. In practical 
feedback systems, imperfections occur such as channel estimation error and feedback delay. These 
inevitable factors degrade the system performance by causing outage lIlTl . |[22l . Therefore, rather than 
using average sum rate as the performance metric, we employ the notion of average goodput Il23l - ll25l 
to incorporate outage probability. Under the multi-cluster subband fading model, we perform analysis on 
the average goodput and the average outage probability with heterogeneous partial feedback. In addition 
to examining the impact of imperfect feedback on multiuser diversity ll26l . ll27l . we also investigate how 
to adapt and optimize the average goodput in the presence of these imperfections. We consider both the 
fixed rate and the variable rate scenarios, and utilize bounding technique and an efficient approximation 
to derive near-optimal strategies. 

To summarize, the contributions of this paper are threefold: a conceptual heterogeneous feedback 
design framework to adapt feedback amount to the underlying channel statistics, a thorough analysis 
of both perfect and imperfect feedback systems under the multi-cluster subband fading model, and 
the development of approximations and near-optimal approaches to adapt and optimize the system 
performance. The rest of the paper is organized as follows. The motivation under the general correlated 
channel model and the development of system model is presented in Section |lll Section |lll] deals with 
perfect feedback, and Section |IV] examines imperfect feedback due to channel estimation error and 
feedback delay. Numerical results are presented in Section |Vl Finally, Section |Vl] concludes the paper. 

II. System Model 

A. Motivation for Heterogeneous Partial Feedback 

This part provides justification for the adaptation of subband size with one cluster of users under the 
general correlated channel model, and motivates the design of heterogeneous partial feedback for the 
multi-cluster scenario in Section III-B I Consider a downlink multiuser OFDMA system with one base 



'An initial treatment of a two-cluster scenario was first presented in 1201 . 
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Station and K users. One cluster of user is assumed in this part and users in this cluster are assumed 
to experience the same frequency selectivity. The system consists of A'^c subcarriers. Hk-n, the frequency 
domain channel transfer function between transmitter and user k at subcarrier n, can be written as: 



L-l 



aiFk^iexp I I , (1) 

1=0 \ c / 

where L is the number of channel taps, for / = 0, . . . , L — 1 represents the channel power delay profile 
and is normalized, i.e., X]£o^ ~ ^' denotes the discrete time channel impulse response, which is 
modeled as complex Gaussian distributed random processes with zero mean and unit variance CJ\f{0, 1) 
and is i.i.d. across k and I. Only fast fading effect is considered in this paper, i.e., the effects of path loss 
and shadowing are assumed to be ideally compensated by power controj^. The received signal of user k 
at subcarrier n can be written as: 

where is the average received power per subcarrier, Sk,n is the transmitted symbol and „ is the 
additive white noise distributed as CA/'(0,cj^ ). From ([T|), it can be shown that H^^ n is distributed as 
CAA(0, 1). The channels at different subcarriers are correlated, and the correlation coefficient between 
subcarriers ni and n2 can be described as follows: 

(TT u \ \^ 2 ( j27r/(n2 -ni)\ 
1=0 ^ ^ 

In general, adjacent subcarriers are highly correlated. In order to reduce feedback needs, subcarriers 
are formed as one resource block, and 77 resource blocks are grouped into one subbancj^. Thus, there are 
N = ^ resource blocks and — subbands ^ . In this manner, each user performs subband-based feedback 

tic f] 

to enable opportunistic scheduling at the transmitter. Since the channels are correlated and there is one 
CQI to represent a given subband, the CQI is a function of the all the individual channels within that 

^This assumption has been employed in 1101 . 1171 , 1271 to simplify the scheduling policy. With the same average SNR, the 
opportunistic scheduling policy is also long-term fair. When different average SNR is assumed, the proportional-fair scheduling 
policy |6| can be utilized. 

'E.g., in LTE, one resource block consists of 12 subcarriers, and one subband can contain 1 to 8 resource blocks 1281 . 
''Throughout the paper, A'c, A'^ and x] are assumed to be a radix 2 number. A more general treatment is possible but this will 
result in edge effects making for more complex notation without much insight. 
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Fig. 1. Comparison of average sum rate for different subband sizes (rj = 1, 2, 4) and partial feedback (M = 2, 4) with respect 
to the number of users. A general correlated channel model is assumed with an exponential power delay profile. (A'c = 256, 
iV = 32, L = 16, (5 = 4, = 10 dB) 



subband. Herein, we employ the following subband (aggregate) average rate S^^r as the functional fo 
IMl, iMl of the CQI for user k at subband r: 

1 '^'^^'^ / PIW |2\ 

' n={r-l)riR, + l ^ "c / 

Each user employs the best-M partial feedback strategy and conveys back the M best CQI values 
selected from Sk^r, 1 < < A detailed description of the best-M strategy can be found in |[T5l . ifTTl . 
|[T9l . After the base station receives feedback, it performs opportunistic scheduling and selects the user 
k for transmission at subband r if user k has the largest CQI at subband r. Also, it is assumed that if no 
user reports CQI for a certain subband, scheduling outage happens and the transmitter does not utilize it 
for transmission. 

Now we demonstrate the need to adapt the subband size to achieve the potential "match" among 
coherence bandwidth, subband size and partial feedback through a simulation example. The channel is 

^This functional form employs the capacity formula and the resulting effective SNR has a geometric mean interpretation. 
Other functional forms of the CQI exist in practical systems such as exponential effective SNR mapping (EESM) f29ll -131 1 and 
mutual information per bit (MMIB) 1321 . 1331 to map the effective SNR to the block-error-rate (BLER) curve. The intuitions 
are similar: to obtain a representative CQI as a single performance measure corresponding to the rate performance. 
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Fig. 2. Illustration of the multi-cluster subband fading channel model for two different clusters with 16 resource blocks. The 
subband sizes equal 2 and 4 for the two different clusters respectively. According to the subband fading model, the channel 
frequency selectivity is flat within each subband, and independent across subbands. The subband sizes can be heterogeneous across 
clusters, and this leads to heterogeneous channel frequency selectivity across clusters. The subband fading model approximates 
the general correlated channel model, and is useful for statistical analysis. 



modeled according to the exponential power delay profile Il36]| - ll38l : of = Yr^[zx75j exp (— |) for 

< I < L, where the parameter 6 is related to the RMS delay spread. The simulation parameters are: 
Nc = 256, N = 32, L = 16, 6 = 4, ^ = 10 dB. The subband size ij can be adjusted and ranges from 

1 to 4 resource blocks. We consider partial feedback with M = 2 and M = 4. The average sum rate of 
the system for different subband sizes and partial feedback with respect to the number of users is shown 
in Fig. [U Under the given coherence bandwidth, several observations can be made. Firstly, the curves 
with 1] = A has the smallest increasing rate because a larger subband size gives a poor representation 
of the channel. Secondly, the curve with r] = 1,M = 2 has the smallest average sum rate because a 
small amount of partial feedback is not enough to reflect the channel quality. Thirdly, the two curves 
r] = 1, M = A and rj = 2, M = 2 possess similar increasing rate. This is because the underlying channel 
is highly correlated within 2 resource blocks and thus having M-best feedback with t] = 2 yields similar 
effect as having 2M-best feedback with rj = 1. From the above observations, ry = 2 matches the frequency 
selectivity and there would be performance loss or waste of feedback resource when a subband size is 
blindly chosen. In a multi-cluster scenario where users in different clusters experience diverse coherence 
bandwidth, this advocates heterogeneous subband size and heterogeneous feedback. 

The general correlated channel model as well as the non-linearity of the CQI, though useful to 
demonstrate the need for heterogeneous feedback, does not lend itself to tractable statistical analysis. 
To develop a tractable analytical framework, an approximated channel model is needed. A widely used 
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model is the block fading model in the frequency domain ||39l . BOl due to its simplicity and capability 
to provide a good approximation to actual physical channels. According to the block fading model, the 
channel frequency selectivity is flat within each block, and independent across blocks 0, iflTl , |[T9l . 
Herein, we generalize the block fading model to the subband fading model for the multi-cluster scenario. 
We assume that users possessing similar frequency selectivity are grouped into a cluster and the subband 
size is perfectly matched to the coherence bandwidth for a given clusteo According to the subband fading 
model, for a given cluster with a perfectly matched subband size, the channel frequency selectivity is flat 
within each subband, and independent across subbands. Fig. |2] demonstrates the subband fading model 
for two different clusters with different subband sizes under a given number of resource blocks. 



B. Multi-Cluster Subband Fading Model 

We now present the multi-cluster subband fading model. Consider a downlink multiuser OFDMA 
system with one base station and G clusters of users. The system consists of resource blocks and the 
total number of users equals K. Users in cluster JCg are indexed by the set /Cg = {1, . . . , /c, . . . , Kg} for 
1 < 5 < G, with \]Cg\ = Kg and ^^^^ Kg = K. In our framework, users in the same cluster group their 
resource blocks into subbands in the same manner while each cluster can potentially employ a different 
grouping which enables the subband size to be heterogeneous between clusters. Denote rjg as the subband 
size for cluster JCg, and r]g G {2*^, 2^, . . . , N}. The ry^'s are ordered such that r/i < • • • < rjc- Based on 
the assumption for rjg, the number of subbands in cluster ICg equals 

Let H^f^, be the frequency domain channel transfer function between transmitter and user k in cluster 
ICg at subband r, where I < k < Kg,l < r < ^. H^^^ is distributed as CA/'(0, 1). According to 
the subband fading model, Hjfj, is assumed to be independent across users and subbands in cluster 
JCg. The feedback for different clusters is at different granularity, and so to model the channel for the 
different clusters of users at the same basic resource block level, some additional notation is needed. Let 
Hjfl = H^^h. with 1 < n < iV denote the resource block based channel transfer function. Then the 
received signals of user k in cluster JCg at resource block n can be represented by: 

4i = ^Hlflsi% + vi% (5) 

*In practical systems, since tiie coiierence bandwidth is determined by the channel statistics which vary on the order of tens 
of seconds or more, the cluster information can be learned and updated through infrequent user feedback. Therefore, the cluster 
is formed dynamically but in a slow way compared to the time variation of the fast fading effect which is on the order of 
milliseconds. 
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where P is the average received power per resource block, s^^^ is the transmitted symbol and f^^^ is 
additive white noise distributed with CA/'(0,cr^). 

Let Zjfl = \hI^I\^ denote the CQI for user k in cluster Kg at subband r. In order to reduce the 
feedback load, users employ the best-M strategy to feed back their CQI. In the basic best-M feedback 
policy, users measure CQI for each resource block at their receiver and feed back the CQI values of the 
M best resource blocks chosen from the total N values. For each resource block, the scheduling policy 
selects the user with the largest CQI among the users who fed back CQI to the transmitter for that resource 
block. However, in our heterogeneous partial feedback framework, since the number of independent CQI 
for cluster /C^ is a fair and reasonable way to allocate the feedback resource is to linearly scale the 
feedback amount for users in cluster ICg. To be specific, user k in ICq (i.e., the cluster with the largest 
subband size) is assumed to feed back the M best CQI selected from {zj^J}, 1 < r < whereas user 
k in ICg conveys the best CQI selected from {zjfl},l < r < ^. After receiving feedback from all 
the clusters, for each resource block the system schedules the user for transmission with the largest CQI. 
It is useful to note that the user feedback is based on the subband level, while the base station schedules 
transmission at the resource block level. 

III. Perfect Feedback 

In this section, the CQI are assumed to be fed back without any errors and the average sum rate is 
employed as the performance metric for system evaluation. We derive a closed form expression for the 
average sum rate in Section IIII-AI for the multi-cluster heterogeneous feedback system. In Section IIII-BI 
we analyze the relationship between the sum rate ratio and the choice of the best-M. 

A. Derivation of Average Sum Rate 

According to the assumption, the CQI zjf^ is i.i.d. across subbands and users, and thus let Fz denote 
the CDF. Because only a subset of the ordered CQI are fed back, from the transmitter's perspective, if it 
receives feedback on a certain resource block from a user, it is likely to be any one of the CQI from the 
ordered subset. We now aim to find the CDF of the CQI seen at the transmitter side as a consequence 
of partial feedback. Let Y^^^ denote the reported CQI viewed at the transmitter for user k in ICg at 
resource block n. Also, let Y^^^ represent the subband-based CQI seen at the transmitter for user k in 
JCg at subband r, then Y^^^ = Y^^)_^-. . The following lemma describes the CDF of Y^^^ (the index n is 
dropped for notational simpUcity), which is denoted by Fy(g) . 
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Lemma 1. The CDF of Y^^^ is given by: 



Fy,.^{x)= UN,M,rj,m){Fz{x)y^ (6) 



m=0 



where the vector = (r/i, • • • ,r]gr ' ' j Vg) and 



-M V « / V"^ 

Proof: The proof is provided in Appendix |A] ■ 
Let /c* demote the selected user at resource block n, then according to the scheduling policy: 

^: = argmax{yi;;,... ,Y^^l... (8) 

where Un = {Kn\- Mn\--- is the set of users who convey feedback for resource block n, 

with \Un^ I = Tg representing the number of users belonging to Un in cluster ICg. It can be easily seen that 
in the full feedback case, i.e., M = Mp = \Un \ = Kg. For the general case when 1 < M < Mp, 
the probability mass function (PMF) of Un is given by: 

nZ^n)= I II C^l I ^<r,<Kg. (9) 




Remark: Only the largest subband size r/c appears in the expression of P(ZYn) instead of the vector rj. 

Zg convey back the ^ 



This is due to our heterogeneous partial feedback design to let users in cluster JCg convey back the 



best CQI out of ^ values. 

Now we turn to determine the conditional CDF of the CQI for the selected user at resource block n, 
conditioned on the set of users providing CQI for that resource block. Since users are equiprobable to 
be scheduled according to the fair scheduling policy, the condition on A;* is not described explicitly, and 
so we denote the conditional CDF as Fx^\u^, where Xn\Un is the conditional CQI of the selected user 
at resource block n. Notice from Lemma [T] that Y^^^ possess a different distribution for different g due 
to the heterogeneous feedback from different clusters. Using order statistics BTI yields i^x„|w„ as: 

G 

Fx„\uS^) = ^{Fyi.^{x)y^. (10) 

9=1 

Then the polynomial form of i*x„|w„ can be obtained, which is stated in the following theorem. 
Theorem 1. The CDF of Fx^\u„ is given by: 

$(M,T7,t) 

Fx^\uSx)= E eG-i(iV,M,r/,T,m)(Fz(x))^"-t^^-'", (11) 

m=0 
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where the vector r ^ (n, • • • , r^, • • • , tg), $(M, r/, r) ^ J2a=i ( - 1 ) , 



Gg(iV,M,r7,r,m) = < 



E Ai(A^,M,T7,r,i)A2(A^,M,r7,T,m-i), 5 = 1 

m 

E Qg-i{N, M, T], T, i)A,+i iN,M,r],T,m-i), 2<g<G 

i=0 



(12) 



Ag{N,M,r],T,m) = < 



(13) 



{Cg{N,M,r],0)y'^, m = 

m^AN,M,r,,0) \{rg + l)£ - m) 

xCg{N, M, T7, i)Ag{N, M,v,T,m-i), I < m < Tg{^M - 1) 
[ {^g{N,M,r^,^M-l)y^, m = Tgi^M-l). 
Proof: The proof is provided in Appendix |A] ■ 
After obtaining the conditional CDF Fx^^\u^, let Cp{M) denote the average sum rate and it can be 
computed using the following procedure. 

1 ^ 

Cp{M) = -^E[log^{l + Xn)] 



(a) 



n=l 

' POO 

Eu / log2(l + px)d(Fx|w(x)) 



V eG-i{N,M,ri,T,m) / log^^l + px)d{Fz{x)) 



(c) 



<S>(M,ri,T) 



G 



N 



J^FiU) QG-iiN,M,rj,T,m)lilp,Y,-^9 

T^O ■m=0 \ 3=1 



m 



(14) 



where p = ^ and F(U) is given by Q. (a) follows from the conditional expectation of Xn\Un and the 
identically distributed property (let X and U represent X„ and Un respectively), (b) follows from ([TTI i in 
Theorem[T] (c) follows from (|9]l, and define Ii{a, h) = log2(l + ax)(i(Fz(x))''. Xi(a, h) is computed 
in Appendix |A] to be: 

6-1 



'+1 



+ 1 



(15) 



where Ei{x) = exjp{—t)t ^dt is the exponential integral function 

The average sum rate for the full feedback is a special case and is given by: 

Cp{Mf)= / log2(l + px)(i(Fz(x))^=Xi(p,if). 
Jo 

Remark: It is noteworthy to mention that the functional form of Cp{M) in ([141) consists of two main 
parts. The first part, which involves ¥{U) and Qg~i{-, •, •, •, •)> accounts for the randomness of the set of 



(16) 
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users who convey feedback as well as the scheduling policy. This part is inherent to the heterogeneous 
partial feedback strategy, and is independent of the system metric for evaluation, such as the average sum 
rate employed in this paper. The second part •) depends on statistical assumption of the underlying 
channel and the system metric, and it is impacted by partial feedback as well. 

B. Sum Rate Ratio and Best-M 

We now examine how to determine the smallest M that results in almost the same performance, in 
terms of average sum rate, as the full feedback case. Applying the same technique as in |[T5l . |fT9l . define 
7p as the spectral efficiency ratio and the problem can be formulated as: 

CpiM*) 

Find the minimum M*, s.t. = — — > 7. (17) 

Cp{Mf) 

The above problem can be numerically solved by substituting the expressions for Cp{M) and Cp{Mf). 
In order to obtain a simpler and tractable relationship between M and K given rj, i.e., the tradeoff 
between the amount of partial feedback and the number of users given existing heterogeneity of channel 
statistics in frequency domain, an approximation is utilized similar to that in llT9l . by observing that 
Zi{a,b) in (ITSl ) is slowly increasing in b with fixed a (This phenomenon is due to the saturation of 
multiuser diversity [43]). Observing X]m=o ''^^ &G-i{N, M,T],T,m) = 1 and employing the binomial 
theorem yields the approximation for the spectral efficiency ratio as: 

7P - 1 - f 1 - ^^^1 • (18) 



N 

From ([TTll and ([TSll . the minimum required M* can be obtained as follows: 



M* > — (1- (1-7)^) . (19) 



VG 

Remark: It can be seen that M* depends on the system parameters (A^, K, 7) as well on the largest 
subband size tjg- It is also a consequence of our heterogeneous partial feedback assumption to let users 
in cluster IC„ convey back the best CQI out of ^ values. This results in the fact that obtaining 

feedback information from users belonging to different clusters have almost the same statistical influence 
on scheduling performance. 

IV. Imperfect Feedback 

After analyzing the heterogeneous partial feedback design with perfect feedback, we turn to examine 
the impact of feedback imperfections in this section. We develop the imperfect feedback model due to 
channel estimation error and feedback delay in Section ITV-AI and investigate the influence of imperfections 
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on two different transmission strategies in Section IIV-B I and IIV-CI Then we propose how to optimize the 
system performance to adapt to the imperfections in Section IIV-DI 

A. Imperfect Feedback Model 

The imperfect feedback model is built upon the subband fading model for the perfect feedback case. To 
differentiate from the notation for the perfect feedback case and focus on the imperfect feedback model, 
the resource block index is dropped. Let denote the frequency domain channel transfer function of 
user k (users in different clusters are not temporally distinguished to avoid notational overload). Due to 
channel estimation error, the user only has its estimated version hk, and the relationship between hk and 
hk can be modeled as: 

hk = hk + Wk, (20) 

where Wk ~ CJ\f{0, a"^^ ) is the channel estimation error. The channel of each resource block is assumed to 
be estimated independently, which yields the channel estimation errors Wk i.i.d. across users and resource 
blocks, i.e., Wk ~ CM{0, a"^). It is clear that the base station makes decision on scheduling and adaptive 
transmission depending on CQI, a function of h^. Thus this information can be outdated due to delay 
between the instant CQI is measured and the actual instant of use for data transmission to the selected 
user. Let be the actual channel transfer function and we employ a first-order Gaussian-Markov model 
||22]| . ||25]| . ||27l to describe the time evolution and to capture the relationship with the delayed version as 
follows: 

hk = ak{hk + Wk) + \Jl- alek, (21) 

where Sk accounts for the innovation noise and is distributed as CJ\f{0, 1). The delay time between 
hk and hk is not explicitly written for notational simplicity, and G [0, 1] is used to model the 
correlation coefficient. Since the feedback delay is mainly caused by the periodic feedback interval and 
processing complexity [25], the innovation noise Sk are i.i.d. across users and a common a is assumed. 
Moreover, Wk and Sk are assumed independent. Therefore, for notational simplicity, the user index k in 
the aforementioned parameters is dropped and Z = is denoted as CQI. 

Let X, X ^nd X represent: the actual CQI of the selected user for transmission, its outdated version, and 
its outdated estimate respectively (x corresponds to X for the perfect feedback case in Section IIII-AI ). 
Notice that the PDF of the outdated estimate x depends on the heterogeneous feedback design and the 
scheduling strategy, whereas the conditional PDF of x|x only depends on a and cr^. Employing the same 
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method in |[26l . |[27l . the conditional PDF is obtained as follows: 

fmi^lx) = ^ exp + I,{alay^), (22) 

where = ■\J^^t^y^^, and /o(-) is the zeroth-order modified Bessel function of the first kind |!42l. 

Since the feedback is imperfect, there are two types of issues that arise. The first is the choice of the 
incorrect user to serve. However, because of the i.i.d nature of the errors this does not compromise the 
fairness and also does not complicate the determination of the CDF. The second problem is that of outage 
because the rate adaptation is made by the base station based on the erroneous CQI. Because of the error 
in the CQI, the rate chosen may exceed the rate that the channel can support and so the base station has 
to take steps to mitigate this effect of outage. A conservative strategy will result in less outage but under 
utilization of the channel while an aggressive strategy will result in good utilization of the channel but 
only for a small fraction of the time. We now present two transmission strategies to address the outage 
issue. 

B. Fix Rate Strategy 

In the fix rate conservative scenario, a system parameter Pq is chosen for rate adaptation, and outage 
results under the following condition: 

Declare outage if : {x < /^olx}- (23) 

The system average goodput is defined as the total average bps/Hz successfully transmitted (23). We 
derive the average goodput and average outage probability for a given choice of system parameter (3q in 
the following procedure. 

Firstly the conditional outage probability is expressed as: 

nx</3o|x) = l-Ql(a«;aVX,«^«V^), (24) 

where Qi(a, b) = Jl^ texp{-^^^)Io{at)dt is the first-order Marcum-Q function EH. Denote i?o(/3o, M) 
as the average goodput for the heterogeneous partial feedback system, which is written according to 
definition: 

Ro{(3o, M) = [E^iw [Hx > /3o|x) log2(l + P/^o)]] ■ (25) 
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Then, from ^ and ([141) . Rq{I3q, M) can be computed as: 

'<l>{M,ri,T) 



y2 @G-i{N,M,r],T,m) Ql{awa^/x,aw\/%)log2{l + pl3o)d{F^{x)y 

Jo 



^8=1 Vn 3 



m=0 

*(M,77,t) / G 



N 



= Y,nU) 0G-i(iV,M,r7,r,m)log2(l + />/3o)X2 /3o,J]-r3-m , (26) 

T^O m=0 y g=l ) 

where X2{a,h) = Qi(a^a-y/x, a^-\/a)(i(F^(x))^. X2{a,h) is computed in Appendix IbI to be: 



(27) 



where = at„a, -d = aw^/a, Q = ^^z^- 

The average outage probabiUty Po(/3oi M) for the heterogeneous partial feedback design can be directly 
calculated from definition and (|26]l as follows: 



Poi(3o,M)=Eu[E^lu[nx<m)]] 

= E^(^) E eG-i(iV,M,r?,T,m) I-X2 /3o,J]-rg-m . (28) 
T^O m=0 y y 3=1 ^9 J J 

The average goodput and average outage probability for the full feedback scenario is a special case 

and is given by: 

Ro{^o,Mf) = log2(l + pPo)MPo, K), 

Po{^o,Mf) = 1-I2{^o,K). (29) 

C. Variable Rate Strategy 

Instead of choosing a conservative system parameter to account for the fix rate scenario as in the 
previous subsection, we consider an approach we refer to as the variable rate strategy. In the variable 
rate scenario, a system parameter /3i is chosen and outage results under the following condition: 

Declare outage if : {x < /^ix\x}, (30) 

where /3i can be regarded as the backoff factor. The system average goodput and average outage 
probability can be derived utilizing the following procedure. 
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Now under the variable rate scenario, the conditional outage probability is expressed as: 



P(x</3ixlx) = l-Qi(a»aV^,a«-VA^)- (31) 

Using the same method as (|25] ) and (|26l ). let Ri{(3i,M) denote the average goodput for the variable rate 
scenario whose expression can be written as follows: 

i?i(/3i, M) = Eu [E^\u inx > /3ix|x) log2(l + pPix)]] 

■f(M,r,,T) I ^ N \ 

= Y,m) 0G-i(A^,M,77,T,m)X3 /3i,J]-r^-m , (32) 

TT^O m=0 \ <7=1 / 

where Z^{a, h) = Qi(a«,a\/i, Ow^/ax) log2(l + pax)d{F ^{x))^ . 
For the full feedback case, the average goodput is given by: 

Ri{pi,MF)=Mlii,K). (33) 

Note that unlike X2{a,h), I^{a,b) can not be written in closed form. Therefore, bounding technique 
and suitable approximation are attractive to find closed form alternatives for X^{a,b). The following 
proposition presents a valid closed form upper bound for T^{a, h) in the low SNR regime. 

Proposition 1. In the low SNR regime, X3(a, h) can be efficiently upper bounded by: 



+ — (— 2i^i(;^,2;2;I^)-±2Fi(l,;^;l;I^) ) ) ), (34) 



where w = Owa, i} = a^-y/a, Q = '^i^l^ , (fi = w'^ + + Q, and 2-^i(-, •; S •) is the Gaussian 
hypergeometric function B2| . 

Proof: The proof is provided in Appendix |B] ■ 
I^^{a, b) is valid and tight especially for the low SNR regime. In order to track ^3(0, b) over the whole 
SNR regimes, we propose the following approximation method by leveraging Jensen's inequality [ 45] . 
Recall the definition of ^3(0, 6) = E[Ql{awa^/x,aw^/(^x)^og2{l + pax)]» where the random variable 
X is defined to have CDF (F^(x))''. Firstly, E[x\ can be computed and is given by: 
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Then plugging (1351 ) into Ql{awa^yx, aw^/ax) log2(l + pax) yields: 



X^(a, b) = Qi (a^a^EM, a^^/^E[^^ log2(l + paE[x]) 



(36) 



Note that T^{a, b) would serve as an upper bound from Jensen's inequaUty if the function of interest 
Ql{al|Ja^/x, a^-^/ox) log2(l + pax) was concave in x. Properties of this function such as monotonicity 
and concavity are of interest and lead to rigorous arguments in support of this bound. If outage does 
not occur, extensive analysis can be carried out due to the well known properties of the log(-) function. 
However, the concavity (or log-concavity) of Qi (a^a-^/x, aw\/Pix) in x (notice that x appears in both 
entries of Qi(-, •)) still remains an important open problem [46 1. Our numerical evidence suggests that 
Ql{awa^/x, aw^/(3lx) log2(l + pl^ix) is concave and monotonically increasing in x for practical choices 
of /?!. For any given /3i preserving the aforementioned property, Jensen's inequality yields an upper bound, 
whose tightness is of interest and discussed in the following proposition. The word practical is used to 
exclude the situation when /3i approaches its maximum 1 which in turn enables Qi(-,-) to dominate 
the goodput to incur extreme outage. This makes intuitive sense according to the definition of average 
goodput. 



Proposition 2. Let {xh} be the family of positive i.i.d. random variables. If Ql{auJa^/x, a^V f^i^) log2(l+ 

p(3ix) is concave and monotonically increasing in x for any given then the Jensen bound is asymp- 
totically tight, i.e., as 6 — )- oo. 



Nonetheless, when the aforementioned property is not preserved (e.g., /3i approaches 1), Jensen's 
inequality does not hold but the expression has been experimentally found to be a good approximation 
and so can still be used. Therefore, (|36] | is denoted as Jensen approximation. We conduct a numerical study 
and demonstrate the tightness of Jensen approximation in Fig. [3l It is observed that the approximation 
method is very tight for moderate (even small) number of users and for all values of Pi G [0, 1], which 
shows its potential in accurately tracking the performance of average goodput. 

Now we calculate the average outage probability. Since it does not involve the log(-) function, it can 





E[Qi(a^a^/X6, a.u,y/ PiXh) log2(l + pPiXh)] 



(37) 



Qi(a^ayEfo], a^y^^Efol) log2(l + p/3iE[x6]) 
Proof: The proof is provided in Appendix |B] 
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Fig. 3. Calculating the average goodput from numerical evaluation {M — 2, 4, 16) and Jensen approximation (M = 16) for 
the variable rate scenario under different p. (TV = 64, r) = (1,4), Ki = K2 = K/2 = 10, a = 0.98, al = 0.01, p = 10 dB, 
and 20 dB) 



be computed into closed forni as follows: 

Pi(/3i,M) = Eu [E^iu [P(X < Pix\x)]] 

HM,V,-r) I ( ^ N 

= Y,m) 0G-i(A^, M,r?,T,m) I I-X4 I A, J] -Tg-m I I , (38) 

where 



Xi{a^h) = I Ql{awa^/x,aw^/ax)d{F^{x))^ 
Jo 

(a) 2b ^^^{-lY j Ql{awax,aw^/ax)exp xdx 



zu = OujO, 1? = a^,^/a, Q = ^pe = + -d"^ + Ce, ^1 = - 'd'^ ^ Q, = - AwH'^. (a) 

follows from change of variables; (b) follows from applying BTl B.48]. 
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In the case of full feedback, the average outage probability Pi{(]i, Mp) becomes: 

Pi(/3i,Afp) = l-X4(/3i,^). (40) 

D. Optimization and Adaptation to Imperfections 

We have obtained the relationship between the system parameter (^o or /3i) and the system average 
goodput, and we now aim to maximize the average goodput by adapting the system parameters. 

Consider the optimization of Ri {f3i , M) to obtain the optimal backoff factor /5J' . It is observed from 
(l32l ) that directly optimizing i?i(/3i, M) is tedious, and a near-optimal method is now proposed to obtain 
I3i . This method is inspired by the results in Section IIII-BI which show that the minimum required M* 
can be chosen to achieve almost the same performance as a system with full feedback. Thus an optimal 
Pi for the full feedback scenario can be optimized first, and then M* is obtained to "match" the system 
performance. Looking again at Fig. |3] with emphasis on different number M of partial feedback, as M 
gets larger, the optimal /3i converges to the full feedback case. In this example, M* = 4 is adequate 
to match the system performance. It is noteworthy to mention that this adaptation philosophy can be 
applied to partial feedback systems wherein system parameters are optimized according to full feedback 
assumption first and minimum required partial feedback is chosen subsequently. 

Note that a closed form approximation has been obtained to track i?i(/3i, Mp) in Section ITV-C[ which 
is denoted as Rf{(3i,Mp) = I^{(3i,K). The following proposition demonstrates the optimal property 
of /3i when optimizing Ri{f3i,Mp). 

Proposition 3. There exists a unique global optimal (3i that maximizes R^{(3i, Mp). 

Proof: The proof is provided in Appendix |B] ■ 
From the above analysis, the optimization strategy can be described as: 

ft* = arg max rHBi, Mp) o:i aig max RUBi.Mp). (41) 

0</3i<l 0</3i<l 

Since it is proved in Proposition |3] that Rf-{l3i, Mp) is quasiconcave ll45l in numerical approach 
such as Newton-Raphson method can be applied to obtain f3^. As discussed before, once (3^ is found, 
the minimum required M* can be obtained by solving dTTl ) or relying on ([T9] l. 

The same strategy can be carried over to the optimization of /3o, which is presented as follows: 

= arg max Ro{Po, Mp). (42) 

The impact of imperfections on system parameter adaptation, and the comparison between the fixed rate 
and variable rate strategies will be examined through simulations in Section |V-B[ 



TO APPEAR IN IEEE TRANSACTIONS ON SIGNAL PROCESSING 



20 



^— Exact minimum IVI (y = 0.99) 
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(b) 



Fig. 4. The required minimum A4 for heterogeneous perfect feedback design: (a) Comparison of the required minimum M 
between numerically solving ( I17t and using approximation ( I19t under different 7 with respect to the number of users; (A'^ = 64, 
77 — (1, 4), Ki — K2 ~ K/2, p = 10 dB) (b) Computing the required minimum AI with respect to different number of users 
when varying the ratio of the number of users in cluster /Ci. (TV = 64, rj — (1,4), p = 10 dB, 7 = 0.99) 



V. Numerical Results 
In this section, we conduct a numerical study to verify the results developed and to draw some insight. 

A. Perfect Feedback Scenario 

The number of resource blocks N is assumed to be 64 for simulations throughout this section. We first 
consider a 2-cluster system. Fig. |4] (a) plots the minimum required M obtained by directly solving ([TT] ) 
and alternatively by the approximation ([T9] l for two thresholds: 7 = 0.99 and 0.9. Note that the result 
from ([T9I ) is rounded with the ceiling function since the required M is an integer The other simulation 
parameters are 77 = (1,4) (i.e., Mp = 16), and p = 10 dB. It is observed that the results from the 
approximate expression matches quite well with the exact computation. The question of whether the 
required M is sensitive to the partition of users in the system is examined in Fig. |4] (b) wherein the ratio 
of the number of users in cluster /Ci is changed and the minimum required M is depicted for different 
total number of users with threshold 7 = 0.99. Interestingly, the result turns out to be "uniform". As 
discussed in Section Jill it is due to the heterogenous feedback design assumption to let users in cluster 
/Ci consume (4M in this simulation) feedback which results in the fact that obtaining feedback 

information from users belonging to different clusters have almost the same influence on scheduling 
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Fig. 5. Comparison of tlie average sum rate for a 4-cluster system under different feedback strategies with respect to the 
number of users (iV = 64, = (1, 2, 4, 8), Ki = K2 = = Ki = M = 2, 4, p = 10 dB) 



perfomiance. Therefore, the representative simulation setup Kg = K/G can be employed when the 
system performance metric is investigated with respect to the total number of users. 

We now consider a 4-cluster system with subband size vector r] = (1,2,4,8) (i.e., Mp = 8). Fig. 
[5] demonstrates the benefit of using heterogeneous feedback design. One of the competing strategies is 
also heterogeneous, but treats users from each cluster separately. In particular, the system firstly clusters 
the users based on their channel statistics, and then serves the clusters one by one requiring feedback 
only from the served cluster of users. In this way, the feedback amount is varying over time depending 
on the partition of users. This strategy is denoted as separate heterogeneous feedback compared to our 
joint heterogeneous feedback design. The other competing strategies are homogeneous without taking 
advantage of the channel statistics of different users. To maintain at least the same feedback amount for fair 
comparison, each user in the homogeneous case is assumed to feed back rE^=i CQI values. Two 

subband sizes are assumed for the homogeneous feedback. It is clear that for the homogeneous case, users 
in cluster JCi have more independent feedback while users in cluster /C4 suffer from redundant feedback. 
The average sum rate for two different values of M are shown in Fig. |5] The separate heterogenous 
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Fig. 6. Comparison of fixed rate and variable rate strategies under normalized parameter /3 (/? = /3i = /3o/10) for different 
number of users K (N = 64, a = 0.98, al = 0.01, p = 10 dB): (a) Comparison of average goodput; (b) Comparison of 
average outage probability. 



feedback is observed to have the worst performance from a sum rate perspective because it does not fully 
exploit multiuser diversity, but it consumes the least feedback. Our joint heterogenous feedback design is 
shown to perform much better than the two homogeneous strategies for the 4-cluster system. It is due to 
the fact that by considering the existing heterogeneity among users, the proposed heterogeneous design 
can make the best use of the degrees of freedom in the frequency domain in order to enhance the system 
performance as well as reduce feedback needs. 

B. Imperfect Feedback Scenario 

Fig. |6] exhibits the comparison between the fix rate and variable rate outage scenarios as well as the 
effect of the number of users on the optimization of Pq and (3i . In order to show the system performance 
of the two scenarios in one figure, a normalized parameter (3 is defined. While examining the variable 
rate plots /3 = /3i, and when considering the fixed rate plots /3 = /3o/10. The system parameters are: 
a = 0.98, (7^ = 0.01, and p = 10 dB. It can be seen that for both scenarios, larger number of users 
K yields better system performance, i.e., higher average goodput and lower average outage probability. 
This is a consequence of increased multiuser diversity gain to combat the imperfections in the feedback 
system. 

Fig. |7] illustrates the effect of channel estimation error cr^ and feedback delay a on the optimal value 
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of /3o and (3i. Here cx^ is varied from to 0.1, and a is varied from 0.9 to 0.99 in steps of 0.005. It 
can be observed from the changing profiles that both the optimal values of /3o and /3i get smaller as the 
imperfections become worse. Therefore, the system should adjust the system parameters to adapt to the 
encountered imperfections. 

Now we consider the adaptation of system parameters (/3o or and partial feedback in a 4-cluster 
heterogeneous feedback system. The system parameters are: r/ = (1,2,4,8), a = 0.98, cr^ = 0.01, and 
p = 10 dB. For both transmission strategies and for a given number of users K, the optimal value of /3q 
or /3* is first optimized according to the full feedback case discussed in Section ITV-D I Then, a minimum 
required M* is obtained by matching the system performance to that in the full feedback case. Fig. [8] 
demonstrates the average goodput for both transmission strategies with M* and (or /3*). We observe 
that there is almost a constant performance gain for the variable rate strategy compared with the fix rate 
one. This is due to the fact that for the variable rate scenario, the system is adapting the transmission 
parameters conditioned on the past memory even if it is the outdated one. If the channel estimation error 
and feedback delay are not severe, the imperfections can be compensated by multiplying with the backoff 
factor and relying on the past feedback. 



TO APPEAR IN IEEE TRANSACTIONS ON SIGNAL PROCESSING 



24 



5 




4 8 12 16 20 24 28 32 36 40 

Number of Users 



Fig. 8. Comparison of the average goodput for a 4-cluster system with fix rate and variable rate strategies using optimized 
and PI. The average goodput is calculated using the best-M partial feedback scheme when the minimum required Af* is 
computed after obtaining /Jq or (A^ = 64, t] = (1, 2, 4, 8), Ki = K2 = K3 = K4. = K/A, a = 0.98, ai = 0.01, p = 10 
dB) 

VI. Conclusion 

In this paper, we propose and analyze a heterogeneous feedback design adapting the feedback resource 
according to users' frequency domain channel statistics. Under the general correlated channel model, we 
demonstrate the gain by achieving the potential match among coherence bandwidth, subband size and 
partial feedback. To facilitate statistical analysis, we employ the subband fading model for the multi- 
cluster heterogeneous feedback system. We derive a closed form expression of the average sum rate 
under perfect partial feedback assumption, and provide a method to obtain the minimum heterogeneous 
partial feedback required to obtain performance comparable to a scheme using full feedback. We also 
analyze the effect of imperfections on the heterogeneous partial feedback system. We obtain a closed 
form expression for the average goodput of the fix rate scenario, and utilize a bounding technique and 
tight approximation to track the performance of the variable rate scenario. Methods adapting the system 
parameters to maximize the average system goodput are proposed. The heterogeneous feedback design 
is shown to outperform the homogeneous one with the same feedback resource. With imperfections, the 
system adjusting the transmission strategy and the amount of partial feedback is shown to yield better 
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performance. The developed analysis provides a theoretical reference to understand the approximate 
behavior of the proposed heterogeneous feedback system and its interplay with practical imperfections. 
Dealing with the general channel correlation and the corresponding nonlinear nature of the CQI are 
interesting directions for the heterogeneous feedback system. 

Appendix A 

Proof (sketch) of Lemma \^ The methodology is an extension of the work in [19] which deals with 



the homogeneous feedback case with one cluster of users and one specific subband size 

tiber of rej 



Let denote the CDF of vj:^^ = Y,^^\ Substituting the subband size — and the number of reported 



CQI for user k in cluster /Co makes K^o satisfy ©. It can be shown that Fv>(s) (x) = f(Y,)^' < 

rig a Tj. ' k 

x) = ^{Y^^rj^-. < x) = Fy[g){x), which concludes the proof. 

Proof of TheoremU} Substituting the expressions of Fvx,) from Lemma [T] and combining (ITOl ) yield: 



Fx^\uA^) = {Fz{x))'^^-^^:^^ll 



{Fzix)y 



3=1 



Applying EH 0.314] to a finite-order power series in (@3), (^X]rn=o ^"(Fzixj)"^'^ j can be expressed 

as X]m=o' ^ ^"^^Fzix))^'^^ ' ^^'^''^ ^^'^ expression for Ag{N, M,ri,T,m) is described in Theorem [T] 
Note that the coefficients of (^p^^^^-^^ can be computed in a recursive manner. 

Then we employ ||48l 0.316] for the multiplication of power series. For g = I, Qi{N, M, r], r, m) can 
be calculated from Ai(A^, M, r], r, m) and A2(A^, M, rj, r, m) as: 

m 

ei{N,M,r],T,m) = ^ Ai(iV, M, 77, r, i)A2(A^, M, 77, r, m - i). 

i=0 

For 2 < g < G, Qg{N, M, rj, r, m) can be computed from Qg^i{N, M, t], t, m) and Kgj^i{N, M, r], r, m) 
in the following manner: 

m 

Qg{N, M, ri, T,m) = Y^ e<,_i(iV, M, rj, r, i)Ag+iiN, M, r},T,m-i). 
This concludes the proof. 

Derivation ofIi{a, h): From the definition of Z, Fz{x) = 1 — exp(— x) and fz{x) = exp(— x). Thus 
d{Fz{x)f = h{Fz{x)f-^fz{x)dx = 6 ^^^^ (^7^) (-1)^ exp(-(£ + where the last equality 

follows from the binomial theorem. 

Therefore, \og^{l + ax)d{Fz{x)f = ^^-^ (^-i) (-if /^^ ln(l + ax) exp(-(£ + l)x)dx. 
Applying |48. 4.337.2] yields 
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Appendix B 

Derivation 07X2(0, b): From the definition of Z, it can be shown that F^{x) = 1 — exp(— ^}^2 x) and 



fzi^) - 1 


1 

— cr^ 




l2ia,b) = 
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-a2)ln2 
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-a2)ln2 



b-l 
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6-1 
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£=0 

b-i 

E 



£=0 



6-1 



6-1 



6-1 



1) / Ql{auJax,aw^/a)ex.p 



{i + l)x^ 



xdx 



(-1)'^ ( Qi(0,^)+exp( 



1 



(-ir- exp(-— ) + exp( 



2(tu2 + Q) 
2{w^ + Ci) 



)(i-Qi(o,^5l^)) 



)(l-exp( 



2(1X72 + 0) 

(44) 

where w = a^a, -i? = a^-^/a, = ^^r^- (a) is obtained by substituting the expression of (i(F^(x))^ 
and using change of variables; (b) follows from applying f4T, B.18]; (c) follows from using the fact that 

Qi(0,^) = exp(-f). 



)) 



Proof of Proposition \J} 



X3(a,6) 



(a) 



< 
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(1 


-a2)ln2 
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-a2)ln2 
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6-1 
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b-i 
e=o 



-1) / Qi(a^a-v/x, au,-v/ax) ln(l + /oax) exp 

(^+ l)x2 



1 - 0-2 



da; 



(— 1) / Qi(au,ax, a^-y/ax) exp 



1 + 



^92/^72 



4^2 

f] 



2F1 



X dx 



1.1. 1^!^! 



4^2 



-2i^l I l,rr;l 



(45) 



where = a^a, = a^-y/a, = ttJ^> = + 'i?^ + 0, and 2-^i(-, •; S •) is the Gaussian 
hypergeometric function P2l . (a) is obtained by substituting the expression of d{F^{x))^; (b) follows 
from the fact that when /) ^ 1, pax is a tight upper bound for ln(l + pax); note that change of variables 
are used; (c) follows from applying P7] B.60]. 

Proof (sketch) of Proposition^ Define s{xb) — Qi(a«;a\/Xfei olu,\/ PiXb) log2(l + pfiiXb)]- Firstly it 
must be shown that g^^^^^y^ converges to 1 in probability. For Ve > 0, it is now shown that: 



s{Xb) 



■mm 



> e 



sixb) - s{nu) 



(a)^^f s{\xb-mb]\) 



> e 



> e 



0, 



(46) 
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where (a) follows from the concave and monotonically increasing property of s(-): \s{x) — s{y)\ < 
s{\x — y\); (b) follows from the asymptotic scaling rate of E[xb\ and \xb — '^[Xb]\, and the utilization of 
the Chebyshev's inequality. From extreme value theory and asymptotic analysis of order statistics [41], 
ll43l . it is known that the tail behavior of Xb converges to type 3 Gumbel distribution, which enables 
^[Xb] to scale as log 6 and \xb — ^[Xb]\ to scale as log log 6. 

Then a method similar to that in Q can be employed to prove the uniformly integrable property P9l of 



si^fi^b ] ) ' combining the above property along with the convergence in probability leads to convergence 



in the mean P9i . which concludes the proof. 

Proof of Proposition^ It must be shown that Z^{l3i,K) = Qi ^Q^a-Y/E[x], aiu-\//3iIE[x]^ log2(l + 
p/3iE[x]) is strictly quasiconcave in j3i. 

This property can be proved by log-concavity [45 1. It is shown in ||46l . |[50l that Qi{\/a, Vb) is log- 
concave in 6 G [0, oo) for a > 0. Also, log(l + b) is concave thus log-concave in 6 G [0,oo). Since 
log-concavity is maintained in multiplication, Qi(\/a, Vb) log(l + b) is log-concave in 5 G [0, oo). From 
the definition of K), it is now proved to be log-concave in /3i G [0, oo) since E[x] is irrelevant 

to /?!. Therefore, it is quasiconcave in /3i G [0, oo) because log-concave functions are also quasiconcave. 

In addition, it is clear that lim li^(8i,K) = 0. Also, it is now shown that: 

0< lim li'{f3i,K) 
< lim exp log2(l + p/3iE[x]) 

S Um ' .exp (j-^vmi]-'..aVmi] = „, (47) 

where (a) follows from the upper bound Qi(a, 5) < exp ^— ^ for 6 > a > (47); (b) follows 
from applying L'Hospital's rule. Therefore, there exists a unique global optimal /3i which maximizes 
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